Quantum lightcone fluctuations in theories with extra dimensions 



Hongwei Yu * and L. H. Ford ' 

Institute of Cosmology, Department of Physics and Astronomy 
Tufts University, Medford, MA 02155, USA 
(February 7, 2008) 

Abstract 

The effects of small extra dimensions upon quantum fluctuations of the light- 
cone are examined. We argue that compactified extra dimensions modify the 
quantum fluctuations of gravitational field so as to induce lightcone fluctua- 
tions. This phenomenon can be viewed as being related to the Casimir effect. 
The observable manifestation of the lightcone fluctuations is broadening of 
spectral lines from distant sources. In this paper, we further develop the for- 
malism used to describe the lightcone fluctuations, and then perform explicit 
calculations for several models with flat extra dimensions. In the case of one 
extra compactified dimension, we find a large effect which places severe con- 
straints on such models. When there is more than one compactified dimension, 
the effect is much weaker and does not place a meaningful constraint. We also 
discuss some brane worlds scenarios, in which gravitons satisfy Dirichlet or 
Neumann boundary conditions on parallel four-dimensional branes, separated 
by one or more flat extra dimensions. 
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I. INTRODUCTION 



One of the most challenging problems in modern physics is the unification of the grav- 
itational interaction with other known interactions in nature. Many attempts, from early 
Kaluza-Klein theory |]J to present supergravity and superstring theories, all involve go- 
ing to higher dimensions and postulating the existence of extra spatial dimensions. The 
modern view regarding these postulated extra dimensions is that they are a physical re- 
ality, rather than merely a technical intermediate device for obtaining rather complicated 
four- dimensional theories from simpler Lagrangians in higher dimensions. If these extra 
dimensions really exist, one must explain why they are not seen. The usual answer is that 
they curl into an extremely small compactified manifold, possibly as small as the Planck 
length scale, l v i = 1.6 x 10~ 33 cm. Therefore low-energy physics should be insensitive to 
them until distances of the compactification scale are being probed. In general, one has the 
possibility of observing the presence of the extra dimensions in a scattering experiment in 
which energies greater than that associated with the compactification scale are achieved. 
Various upper bounds have been put on the size of possible extra dimensions F° r 
example, an upper bound of ~ 1 Tev was given in orbifold compactifications of superstrings 
Q. However, if only gravity propagates in the extra dimensions, the upper bound can be 
much larger. A recent proposal is that the fundamental scale of quantum gravity can be as 
low as few Tev and the observed weakness of gravity is the result of large extra dimensions 
in which only gravity can propagate |J. This scenario could be realized in the context of 
several string models || in which one has a set of three-branes (3+1 dimensional spacetime) 
in the entire spacetime with extra dimensions. The Standard Model particles are confined 
to one of the branes, while gravitons propagate freely in the entire bulk. The size of extra 
dimensions could then be as large as 1 mm in this type of model. Extra dimensions of 
sufficiently large size may manifest themselves in particle colliders and in the possible 
deviation from Newton's law at short distances ||, and they may also have implications in 
gauge unification |J and cosmology [fjlj. 

However, a question arises naturally as to whether there are any lower bounds on the 
sizes of extra dimensions. It is the common belief that the existence of extra dimensions 
has no effect on low-energy physics as long as they are extremely small. Recently, we have 
argued that this is not the case, because of lightcone fluctuations arising from the quantum 
gravitational vacuum fluctuations due to compactification of spatial dimensions [11]; HI- ^n 
explicit calculation was carried out in the five-dimensional prototypical Kaluza-Klein model 
which showed that the periodic compactification of the extra spatial dimensions gives rise 
to stochastic fluctuations in the apparent speed of light which grow as the compactifica- 
tion scale decreases and are in principle observable. Basically, the smaller the size of the 
compactified dimensions , the larger are the fluctuations that result. This is closely related 
to the Casimir effect, the vacuum energy occurring whenever boundary conditions are im- 
posed on a quantum field. The gravitational Casimir energy in the five-dimensional case 
with one compactified spatial dimension was studied in , where a nonzero energy density 
was found, which tends to make the extra dimension contract. This raises the question of 
stability of the extra dimensions. It is possible, however, that the Casimir energy arising 
from the quantum gravitational field and other matter fields may be made to cancel each 
other |14[], thus stabilizing the extra dimensions. Quantum lightcone fluctuations due to the 
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compactification of spatial dimensions |IT,|I2,15|, although similar in nature to the Casimir 



effect, come solely from gravitons. Hence, no similar cancellation is to be expected. 

In an earlier work fll^fl , we examined lightcone fluctuation in a five- dimensional model 
with periodic compactification. We found that there seem to be observable effects which 
essentially rule out this model. That is, we derived a lower bound on the compactification 
scale which is larger than upper bounds derived from other considerations. The purpose of 
the present paper is threefold: to further develop the basic formalism, to provide more details 
of the five- dimensional calculation, and to extend the analysis to some higher dimensional 
models. 

In Section II, we give a brief review of the formalism, discuss the observability of lightcone 
fluctuations, and derive the graviton two-point functions in the transverse tracefree gauge 
for spacetimes with an arbitrary number of dimensions (detailed calculations will given 
in the Appendix). In Section III, we examine light cone fluctuations in spacetimes with 
extra dimensions periodically compactified into a torus. The five- dimensional prototypical 
Kaluza-Klein model, of which some results have already been reported, will be studied in 
great detail. Higher dimensions up through 11 will also be discussed. Section IV deals with 
light cone fluctuations in the brane-world scenario, as motivated by a recent proposal of 
extra dimensions of macroscopic size. We will summarize and conclude in Section V. 



II. OBSERVABILITY OF LIGHT CONE FLUCTUATIONS AND THE 
GRAVITON TWO-POINT FUNCTION IN THE TT GAUGE 



To begin, let us examine a d = 4 + n dimensional spacetime with n extra dimensions. 
Consider a flat background spacetime with a linearized perturbation h^ v propagating upon 
it , so the spacetime metric may be written as ds 2 = (rj^ + h^dx^dx" = dt 2 — c/x 2 + 
h^dx^dx" , where the indices /z, v run through 0, 1, 2, 3, 3 + n. Let <r(x, x') be one half of 
the squared geodesic distance between a pair of spacetime points x and x', and o"o(x, x') be 
the corresponding quantity in the flat background . In the presence of a linearized metric 
perturbation, h^ u , we may expand o = <r + + Oih 2 ^) . Here o\ is first order in h^ u . If 
we quantize h^, then quantum gravitational vacuum fluctuations will lead to fluctuations 
in the geodesic separation, and therefore induce lightcone fluctuations. In particular, we 
have (a 2 ) ^ 0, since o\ becomes a quantum operator when the metric perturbations are 
quantized. The quantum lightcone fluctuations give rise to stochastic fluctuations in the 
speed of light, which may produce an observable time delay or advance At in the arrival 
times of pulses. 



A. Observability of light cone fluctuation 



Here, we shall discuss how light cone fluctuations characterized by (a 2 ) are related to 
physical observable quantities. For this purpose, let us consider the propagation of light 
pulses between a source and a detector separated by a distance r on a flat background with 
quantized linear perturbations. For a pulse which is delayed or advanced by time At, which 
is much less than r, one finds 
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a = a + a 1 + .... = -[(r + At) 2 - r 2 ] « rAt . (1) 

Square the above equation and take the average over a given quantum state of gravitons \(j>) 
(e.g. the vacuum states associated with compactification of spatial dimensions ), 

A^IMW (2) 

This result is, however, divergent due to the formal divergence of (<ft\<jf\<f)) . One can define 
an observable At f, s by subtracting from Eq. (0) the corresponding quantity, At 2 ,, for the 
vacuum state as follows 

44Hil ;_ 4 ^»H*ii 5 MW. (3 ) 

Here we take the absolute value of the difference between At 2 and At 2 ,, because the observ- 
able quantity At 2 fes has to be a positive real number. Note that we can also get this result 
from the gravitational quantum average of the retarded Green's function (G ret (x,x')) when 
(&i)r > ||15|| . Therefore, the root-mean-squared deviation from the classical propagation 



time is given by 



M<*, = iMM. . (4) 



At this point, a question may arise as to whether the formal procedure of taking the 
absolute value in deriving the relation between At and (a 2 )/?, Eq. (f§), is a reasonable one, or 
whether a meaningful relation between At and {(t\)r can be established only when (cx 2 )^ > 0. 
We shall argue that a result essentially the same as Eq. @j can be obtained by the following 
analysis which avoids the sign problem. Instead of squaring the Eq. (|l|), we take the fourth 
power of both sides and average over a quantum vacuum state to yield 

At <P = r 4 • ( 5 ) 

We can regularize {4>\af\4>) by normal ordering and define 

{MWn = (01 : of : |0) • (6) 
For a free field ip and a quantum vacuum state one can show by use of Wick's theorem that 

<:V 4 :>=3<:^ 2 :) 2 =3<^>|. (7) 

Therefore one finds that 

r r 

This is essentially the same as Eq. (Q), apart from a dimensionless factor of order unity. We 
should note that there may be other ways to define the quartic operator, af. One possibility 
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is to let a\ = (: o[ :) , provided that the integrals involved can be defined. Both of these 



definitions were discussed in Ref. ||16|| . There it was found in some model cases that the two 
definitions yield the same result, apart from numerical factors of order unity, which are not 
important for the present purposes. 

Note that At is the ensemble averaged deviation, not necessarily the expected variation 
in flight time, 5t, of two pulses emitted close together in time. The latter is given by At only 
when the correlation time between successive pulses is less than the time separation of the 
pulses. This can be understood physically as due to the fact that the gravitational field may 
not fluctuate significantly in the interval between the two pulses. This point is discussed 
in detail in Ref. |I7|| . These stochastic fluctuations in the apparent velocity of light arising 
from quantum gravitational fluctuations are in principle observable, since they may lead to 
a spread in the arrival times of pulses from distant astrophysical sources, or the broadening 
of the spectral lines. Lightcone fluctuations and their possible astrophysical observability 
have been recently discussed in a somewhat different framework in Refs. [|T3,|T^]. 



B. An Alternative Derivation of At 



In order to find At in a particular situation, we need to calculate the quantum expectation 
value (u1)r in any chosen quantum state \ip), which can be shown to be given by [] flT5| , |TTf 



{al) R = \{Arf C dr P dr' nVn"n" G* vpa (x, x') . (9) 

O Jr Jr 



Here dr = |dx|, Ar = r% — tq and = dx^ /dr. The integration is taken along the null 
geodesic connecting two points x and x', and 

G %pa{ x , x ') = (^\h^(x)h p(T (x') + h^(x')h pa (x)\^) (10) 

is the graviton Hadamard function, understood to be suitably renormalized. The gauge 
invariance of At, as given by Eq. (|]), has been analyzed recently |p"T| . 

In this subsection, we wish to rederive At using the geodesic deviation equation. This 
derivation allows us to see the gauge invariance more clearly, and to discuss the issue of 
Lorentz invariance of lightcone fluctuations. Let us consider a pair of timelike geodesies 
with tangent vector u^, and n M as a unit spacelike vector pointing from one geodesic to the 
other (See Fig. 1). 



Although the derivations there were given in 3+1 dimensions, the generalization to arbitrary 
dimensions is straightforward. 
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FIGURES 




FIG. 1. A pair of nearby timelike geodesies. Here is a tangent vector along the geodesic, 
while is a unit spacelike vector pointing from one geodesic to the other. 



The geodesic deviation equation is given by 



dr 2 



(11) 



where R^ u g is the Riemann tensor. The relative acceleration per unit proper length of 
particles on the neighboring geodesies is 

D 2 n^ 

a = n^-^- = -R mv ^u a n v vP . (12) 

Thus if ds is the spatial distance between the two particles, then a ds is their relative 
acceleration. It follows that the relative change in displacement of the two particles after a 
proper time T is 



ds dr dr'a(r',0), (13) 
J o Jo 

Now consider the case of two observers (particles) separated by a finite initial distance s as 
illustrated in Fig. (0). 



ds ■ ds 



So 
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FIG. 2. Two timelike geodesies separated by a finite interval containing an infinite number of 
nearby geodesies . 

We can find the relative change in displacement of these two observers by integrating on 

s: 

As = f ° ds ( T dr f dr'air', 0) . (14) 
Jo Jo Jo 

This is the relative displacement measured at the same moment of proper time for both 
observers. 

Let us now consider a light signal sent from one observer to the other. If a = 0, the 
distance traveled by the light ray is so- When a^O, this distance becomes so + As, where 
now 

rso rs rr 

As = / ds dr dr'afr'^s) . (15) 
Jo Jo Jo 

Here the under-braced integral is the displacement per unit s of a pair of observers at 
a distance s from the source. The domain of the final two integrations is illustrated in 
Fig- ©. 

T 1 




FIG. 3. The displacement, As, between a source and a detector is given by an integration 
within the triangular region. 

If gravity is quantized, the Riemann tensor will fluctuate around an average value of 
zero due to quantum gravitational vacuum fluctuations. This leads to (a) = 0, and hence 
(As) = 0. Notice here that a becomes a quantum operator when metric perturbations are 
quantized. However, in general, ((As) 2 ) ^ 0, and we have 

rso rso rsi rrx rS2 PT2 

((As) 2 )=/ d Sl / ds 2 / dn dr[ dr 2 dr^(a{r' 1 ,s 1 )a(T^s 2 )). (16) 
Jo Jo Jo Jo Jo Jo 

Thus the root-mean-squared fluctuation in the flight path is \J ((As) 2 ), which can also be 
understood as a fluctuation in the speed of light. It entails an intrinsic quantum uncertainty 
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in the measurement of distance. Therefore, spacetime becomes fuzzy at a scale character- 
ized by \J ((As) 2 ). The integrand in Eq. ([16]) is obviously invariant under any coordinate 
transformation while the integral is gauge invariant within the linear approximation. 

We now wish to show that this gauge-invariant quantity is the same as Eq. (|4]) when 
calculated in the transverse-tracefree (TT) gauge. Choose a coordinate system where the 
source and the detector are both at rest, and suppose that the light ray propagates in the 
x-direction, then we have 

u" = (1,0,0,0), (17) 
n" = (0,l,0,0), (18) 



and 



ot — Rxtxt — —-h xx ,tt • (19) 
Substitution of the above results into Eq. (|16|) leads to 

/f pf /"^l / 1 "2'2 P^2 

dx l dx 2 dti dt[ dt 2 dt' 2 (a{t' 1 ,s 1 )a{t' 2 ,s 2 )) 
Jo Jo Jo Jo Jo 

1 r r r T 1 
= - / dx x \ dx 2 (h xx (x 1 ,xi)h xx (x 2 ,x 2 )) = -(a 2 ) , (20) 
4 jo Jo r 

where we have set Sq = r and used the fact that along the light ray x = t. Thus, one has 



v (*1> / 

At = ^ = V<(A*) 2 > (21) 

which also demonstrates the gauge-invariance of At. 

Now we wish to discuss the rather subtle issue of the relation of lightcone fluctuations 
to Lorentz symmetry. It is sometimes argued that lightcone fluctuations are incompatible 
with Lorentz invariance. The most dramatic illustration of this arises when a time advance 
occurs, that is, when a pulse propagates outside of the classical lightcone. In a Lorentz 
invariant theory, there will exist a frame of reference in which the causal order of emission 
and detection is inverted, so the pulse is seen to be detected before it was emitted. Thus 
the lightcone fluctuation phenomenon, if it is to exist at all, seems to be incompatible with 
strict Lorentz invariance. 

Our view of the situation is the following: lightcone fluctuations respect Lorentz sym- 
metry on the average, but not in individual measurements. The symmetry on the average 
insures that the mean lightcone be that of classical Minkowski spacetime. The average met- 
ric is that of Minkowski spacetime provided that (h^) = 0. However, a particular pulse 
effectively measures a spacetime geometry which is not Minkowskian and not Lorentz in- 
variant. A simple model may help to illustrate this point. Consider a quantum geometry 
consisting of an ensemble of classical Schwarzschild spacetimes, but with both positive and 
negative values for the mass parameter M. (The fact that the M < Schwarzschild space- 
time has a naked singularity at r = need not concern us. For the purpose of this model, 
we can confine our discussion to a region where r 3> \M\.) Suppose that this ensemble has 
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(M) = 0, but (M 2 ) ^ 0. It is well known that light propagation in a M > Schwarzschild 
spacetime can exhibit a time delay relative to what would be expected in flat spacetime. 
This is the basis for the time delay tests of general relativity using radar signals sent near 
the limb of the sun. In the present model, however, the time difference is equally likely to be 
a time advance rather than a time delay. A measurement of the time difference amounts to 
a measurement of M. This model is Lorentz invariant on the average because (M) = and 
the average spacetime is Minkowskian. However, a specific measurement selects a particular 
member of the ensemble, which is generally not Lorentz invariant. 

In addition to the fact that the mean metric is Minkowskian, there is another sense 
in which lightcone fluctuations due to compactification exhibit average Lorentz invari- 
ance. Note that As, and hence At, depends on the Riemann tensor correlation function 
(Rxtxt(xi)Rxtxt(%2)), which is invariant under Lorentz boosts along the x-axis. Thus if we 
were to repeat the above calculations of As in a second frame moving with respect to the 
first, the result will be the same. In both cases one is assuming that the detector is at 
rest relative to the source. This is a reflection of the Lorentz invariance of the spectrum of 
fluctuations, which is exhibited by the compactified flat spacetimes studied in this paper, 
but not by the Schwarzschild spacetime with a fluctuating mass. 

C. Graviton two-point functions 

We shall use a quantization of the linearized gravitational perturbations h^ v which retains 
only physical degrees of freedom. That is, we are going to work in the TT gauge in which 
the gravitational perturbations have only spatial components h^, satisfying the transverse, 
d l hij = 0, and tracefree, h\ = conditions. Here i, j run from lto3 + n = d— 1. These 2d 
conditions remove all of the gauge degrees of freedom and leave \{d 2 — 3d) physical degrees 
of freedom. We have 

hij = E k,Aey(k, A)/ k + H.c.}. (22) 

k,A 

Here H.c. denotes the Hermitian conjugate, A labels the ^(d 2 — 3d) independent polarization 
states, /k is the mode function, and the e /ii/ (k, A) are polarization tensors. (Units in which 
327iGd = 1, where is Newton's constant in d dimensions and in which h = c = 1 will be 
used in this paper.) 

Let us now calculate the Hadamard function, Gfj, up(T (x, x'), for gravitons in the Minkowski 
vacuum state in the transverse tracefree gauge. It follows that 

G m (x, x>) = -^- x J ^ E A)e„(k, AK k 'C*-* O c -M*-0 . ( 23 ) 
The summation of polarization tensors in the transverse tracefree gauge is (See the Appendix 



in Ref. 11 . 



The tensorial argument given there applies in any number of dimensions 
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ey(k, A) e k i(k, A) = 5 ik Sji + 5 u S jk - 5^5^ + kikjk k ki 

A 

-\-kikjd k i + k k ki5ij kik[5j k kik k 5ji kjki5i k kjk k 5u , 



where hi = y . We find that 

Gij k i = 2Fij5ki + 2FkiSij — 2F ik 5ji — 2Fu5j k — 2Fji5i k — 2Fj k Su 
+2H m + 2D(x, x')(Si k Sji + 5 u 5 jk - S^S^) . 

Here D(x,x'), Fij(x,x') and Hij k i(x,x') are functions which are defined as follows 

n»f r T '\ - Re f rf3+nk r ik-(x-x') -Mt-t') 
[ ' ] ~ (27r)3+™y 2cu 



F!}(x,x" 



Re a f)> f rf3+nk r »k-(x-x') -Mt-t') 
(2n)z+" 1 j J 2^ 



and 



^ q a/ / <f * ik-(x-x') -iu;(t-t') 

2cu 5 



These functions are calculated in the Appendix. Let 

J R=|x-x'|, At = t-t'. 
For n = 2m + 1, an odd number of extra dimensions, the results are 

' fori? 2 < At 2 , 



ij 



2(2tt) 



(2m)!! 
i? 2m 



5(o) E 



(2fc+ 1)!! 



R 



2k 



k=0 



(2k)\\(2k + l)(2k - 1) (i? 2 - At 2 ) k 



and 



2(27r) m+2 J I/? 2 " 1-2 

^ 2 (m-j-l)(2j + l)!! 

4 ™C? + l)(2j)!!(2j + l)(2j - 1) (it> 2 - At 2 )^' 



5(0) 



In particular, for n — 1 (d = 5), we have 



Here 
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and 



f VR 2 -At 2 f n2 v. A f 2 

s ®=\^ (34) 



m = (\-^ m - (35) 



For n = 2m, an even number of extra dimensions, we have 



(2vr) m + 2 {R 2 - At 2 )"^ 1 ' 1 ; 



and 



ij 



1 g,g,f (2m-l)!! 
(27r) m + 2 1 H i? 2m 



At + At\ 2 
4R At, 



_ 1 ™ 2 fc - 2 r(A; 
_ ~D2 2^ 



1} 



2A- 



i? 2 



fc=2 



(2fc-l)!! (i^-At 2 )^ 1 



m > 1 



Tj2m 
n ijkl 



_J_ „,f(2m-l)!! 



At 2 
o 7 ^ 



m-2 



E 

fc=3 



J^2m~A 
1 



At /At 2 \ /R + At\ 2 



2AR 3 V i? 2 
' 1 At 



18i? 2 



=-bl 



(2A; — l)(2fc — 3) V^ 2 AR3 Vi?-At 



i? + At\ 2 



;m-j-i)2^ 2 r(j-i) 



i=2 



(2m-l)(2j + l)H (R 2 -At 2 y~ 1 



m>2. 



For the case n = 2 (d = 6), we have for if 



(37) 



(38) 



H, 



1 



(27T) 



MM 
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At 2 At /At 2 



6# 2 8RK3R 2 



1 In 



g + At \ 2 
R-At 



(39) 



The case of n = (d = 4) was given in Ref . fl 1 j 



III. PERIODIC COMPACTIFICATION SCENARIO 



Let us now suppose that the extra n dimensions zi, z n are compactified with periodicity 
lengths Li, L n , namely spatial points Zi and Zi + Lj are identified. For simplicity, we shall 
assume in this paper that L\ — ... — L n — L. The effect of imposition of the periodic 
boundary conditions on the extra dimensions is to restrict the field modes to a discrete set 

/ k = (2uo(2ufL n )-^e l ^- ult) , (40) 

with 
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27P/77, ■ 

k l = — ^, i = l,...,n, m, = 0,±1,±2,±3,.... (41) 

Let us denote the associated vacuum state by \0l)- in order to calculate the gravitational 
vacuum fluctuations due to compactification of extra dimensions, we need the renormalized 
graviton Hadamard function with respect to the vacuum state |0_c), (x,x'), which is 

given by a multiple image sum of the corresponding Hadamard function for the Minkowski f] 
vacuum, G^ up(T : 

n +00 ' 

G *vpa(t, Z h t', Z[) = J[ G ^pa(t, Zi, t', Z\ + 771 jL) . (42) 

i=l rrii=—oo 

Here the prime on the summation indicates that the m; = term is excluded and the 
notation 

(t,x,zx, ..,z n ) = (t,Zi) (43) 

has been adopted. 

We are mainly concerned about how lightcone fluctuations arise in the usual uncom- 
pactified space as a result of compactification of extra dimensions. So we shall examine the 
case of a light ray propagating in one of the uncompactified dimensions. Take the direction 
to be along the x-axis in our four dimensional world, then the relevant graviton two-point 
function is G xxxx , which can be expressed as 



G xxxx {t } x } Zi } t , x , Zj^) 2 



D(t, x } Zi } t , x , z^j 1F xx {t } x, Zi } t , x , Zj^) 



(44) 



~^~H xxxx {t } X, Zii t , X , z^) 

Assuming that the propagation goes from point (a, 0, 0) to point (6, 0, 0), we have 
(° 2 i) = l(b-a) 2 f dx [ b dx'G* xxx (t,x,0, t',x',0,), 

J a J a 

1 nb rb n +°° ' 

= -(b-a) 2 dx dx' Yl G xxxx(t, x, 0, t', x', 0, 0, m x L, rmL) . 

° Ja Ja i=lmi=-oo 

(45) 

With these results, we can in principle calculate lightcone fluctuations in spacetimes with 
arbitrary number of flat extra dimensions. In what follows, we first examine some particular 
cases, then make some observations for the general case. 



By Minkowski we mean flat spacetime with all dimensions uncompactified 
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A. The five dimensional Kaluza-Klein model 



One of the most intriguing and elegant ways of unifying gauge field theories with grav- 
itation is the higher-dimensional generalizations of Kaluza-Klein theory. The original sug- 
gestion of Kaluza and Klein Jl] was that electromagnetism and general relativity could be 
unified by starting with a five- dimensional version of the latter and then somehow arranging 
for the fifth dimension to become unobservable. This idea was further generalized to higher 
dimensions in attempts to unify non-Abelian gauge fields with gravitation and has been ex- 
tensively studied in recent years in the context of supergravity and superstring theories. In 
the course of investigation of new features arising from the introduction of extra dimensions, 
the five-dimensional Kaluza-Klein theory has always been taken as a prototypical model to 
carry out explicit calculations to obtain a basic understanding. This is also our strategy 



here. In this subsection, we will derive the results which were summarized in Ref. [12 



1. Calculation of At 



To begin, let us examine the influence of the compactification of the fifth (extra) di- 
mension on the light propagation in our four dimensional world, by considering a light ray 
traveling along the x-direction from point a to point b, which is perpendicular to the direction 
of compactification. Define 



p = x — x 



b — a 



(46) 



and note the fact that the integration in Eq. (^) is to be carried out along the classical null 
geodesic on which t — t' = p. Then we obtain, after performing the differentiation in both 



D and F xx , 



G xxxx (t, x, 0, 0, 0, t', x, 0, 0, mL') 



4vr 2 
-+ 



7p 4 



(p 2 + m 2 L 2 ) 3 \mL\ 3 (p 2 + m 2 L 2 ) 3 \mL\ 



9p 2 \mL\ 



\mL\ 



(p 2 + m 2 L 2 ) 3 (p 2 + mL 2 f 



(47) 



Thus, we have 



= \ r2 fdx f dx' G* xxx {t,x,0, f,a/,0,), 

J a J a 

1 rb rb +°° ' 

= -r 2 / dx dx' G xxxx (t,x,0, t',a/,0, 0,mL) 

O J a J a ___ , , 



m=— oo 



32n 2 L 



E 

m=l . 



— In 1 + — T - 7 2 , 2^ 

m m z mr y~) A + m z )m 



(48) 



where we have introduced a dimensionless parameter 7 = r/L. We are interested here in 
the case in which 7 > 1. It then follows that the summation is dominated, to the leading 
order, by the second term, 
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where £(3) is the Riemann-zeta function. So, the mean deviation from the classical propa- 
gation time due to the lightcone fluctuations is 



Here we have used the fact that G5 = G^L, and t p i « 5.39 x 10~ 44 s is the Planck time. 

This result reveals that here the mean deviation in the arrival time increases linearly^ 
with r and inversely with the size of the extra dimension. This effect is counter-intuitive in 
the sense that it grows as the size of the compactified dimension decreases. When r is of 
cosmological dimensions and L is sufficiently small, the effect is potentially observable. 

2. Choice of Renormalization and Changes in L. 

Note that here (<j\)r has been renormalized to zero as L — > 00. This is the most natural 
choice of renormalization, corresponding to the effect of the graviton fluctuations vanishing 
in the limit of noncompactified spacetime. This is analogous to setting a Casimir energy 
density to zero in the limit of infinite plate separation, However, if instead of renormalizing 
(of) against the vacuum with respect to L —>■ 00, we take the manifold with compactified 
extra dimensions at some fixed sizes L to have (of )# — 0, then the lightcone fluctuations 
could seem to be renormalized away. The latter renormalization scheme is a logical possibility 
that we can consider, although it is unnatural as there seems to be nothing in the theory 
which picks out a particular finite value of L. In any case, if L is somehow allowed to vary, 
for example, as the universe evolves, then lightcone fluctuations would produce noticeable 
effects, as one could at most set (of)/? = at one point along the path of a light ray. 
It is particularly so, when we try to detect the spread in the arrival times of pulses from 
distant astrophysical sources, where we are looking back in time. Hence significant lightcone 
fluctuations may arise no matter what renormalization scheme one chooses if the size L is 
allowed to vary. 

To get an understanding for the case of a changing L, let us assume that L changes with 
time at an extremely small rate, which is in fact required by experimental data on the time 
evolution of fundamental constants (see, for example, Refs. [PtJ-pSfl). Then the evolution of 
L can be reasonably well approximated by a linear function of time: f\ 



4 In the usual four dimensional case with one compactified spatial dimension At grows linearly 
with the square root of r (see Ref. Jll[|). 

5 Strictly speaking, the functional dependence of L on time should be given by a yet-unknown 
underlying dynamical theory which governs how the extra dimensions evolve. However, the as- 
sumption of a linear dependence is good enough for our purpose of getting a basic idea about how 
the variation of L over time would affect our results. 
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L = Li + at 



(51) 



Using this expression for L and redoing the calculations [], one finds, 



32n 2 L 
C(3)rV 

167T 2 L f 



f m =l ^ L i 



8 1 n A 27 

— In 1 + — - 



(7 2 + m?)m 



(52) 



where L{ is the initial compactification size when the light ray is emitted, Lf = Li + ar is 
the final size at the time of reception and 7 = r/Lf. Here (<j 2 )r is renormalized with respect 
to L — > 00. Therefore, one has for the mean time deviation from the classical propagation 
time 



At 



r 



(53) 



Another possibility, as we have mentioned earlier, is to renormalize (erf) against that 
corresponding to the current size Lf, which implements the idea of setting the renormalized 
(<j 2 )r to be zero if L is fixed always or at least during the propagation of the light. This is 
accomplished by taking the difference between Eq. ( |5^ ) and Eq. (|33) with L being replaced 
by Lf to obtain 



C(3)rV 
16n 2 L 



(54) 



which leads to 



At 



>i 2 ), 



Li 



Li 



■pi 



(55) 



Equations ( p3|) and flog ) demonstrate clearly that no matter what renormalization scheme 
is employed, one gets a nonzero lightcone fluctuation effect as long as L is changing. We 
want to point out here again that renormalizing (of) with respect to L — > 00 is far more 
natural than to a particular finite size Lf, since the latter seems to pick out a preferred size 
Lf without any convincing reason to do so. 



3. Correlation of Pulses 

The fluctuation in the flight time of pulses, At, can apply to the successive wave crests 
of a plane wave. This leads to a broadening of spectral lines from a distant source. Note, 



6 It is worth pointing out here that the graviton two-point function obtained simply by replacing 
the constant L with a changing one is not the exact two-point function that satisfies the appropriate 
equations, but it is a very good leading order approximation provided that a <C 1 
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however, that At is the expected variation in the arrival times of two successive crests only 
when the successive pulses are uncorrelated [17]. To determine the correlation, we need to 
compare |(cr 2 )| and | (cicri) ] - The latter quantity is defined by 



(a 1 a[) = l(Ar) 2 P dr P dr' 

O Jro J ro 



n»n v n p rf Gl„(x,x') 



fiupcr \ 



(56) 



where the r-integration is taken along the mean path of the first pulse, and the r'-integration 
is taken along that of the second pulse. Here we will assume that At <C Ar, so the slopes, 
v and v', of the two mean paths are approximately unity. Thus the two-point function in 
Eq. ( |56|) will be assumed to be evaluated at p = |x— x'| = \r— r'\ and r = \t — t'\ = \r—r'—t \. 
If | (01 o{)| <C \(crf)\, two pulses are uncorrelated, and otherwise they are correlated. 

Let us now suppose the time separation of two pulses is T, and note that the relevant 
graviton two-point function can be expressed as 



G xxxx (t, x, 0, 0, 0, t', x', 0, 0, nL) \ 



4vr 2 



t-t'=p-T 



+ 



(3 3 (3(p 2 + n 2 L 2 ) 2 (3(p 2 + n 2 L 2 ) 
16/V 4/3 2n 2 L 2 



(p 2 + n 2 L 2 f (p 2 + n 2 L 2 ) 2 f3 3 (p 2 + n 2 L 2 ) 



4vr 2 /? 3 4tt 2 ' 



(57) 



where, 



(3(n,p) 

Utilizing the following integration relation 



n 2 L 2 + 2 P T-T 2 ) 1/2 



(58) 



dx I dx'f(x — x') 



(r-p)[f(p) + f(-p)]dp, 



(59) 



one finds that 



(a 1 a[} = A + B, 



(60) 



where 



A = -TT^T. / d P( r -p) 
r 



n=l 

2 oo t 4 

y — 

167r 2 L 3 A/ T 2 



1 



+ 



2\ n 2 - 



J>2 

L 2 



2rT - T 2 



'n 2 + 



L 2 



2rT + T 2 



' rr 



L 2 



(61) 



and 
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2 oo 

B = T^T, I dp(r-p)[f(p,n) + f(-p,n)] 
16vH ^ Jo 



E2 



2 v / n 2 L 2 - T 2 /3(r, n)+P(-r, n) 2 /nL + ^n 2 L 2 -T 2 
tfL 2 + {n 2 L 2 + r 2 ) + nL \nL - y/nHJ^f^ 



1 j /n 2 L 2 + rT + nL(3(r,ri)\ 1 ^ /n 2 L 2 - rT + nL(3(-r,n) 



nL \n 2 L 2 + rT — nLfi(r,n) J nL \n 2 L 2 — rT — nL(3(—r,n 



(62) 



A few things are to be noticed here: (1) We need to drop those terms in A and B when 
the square root is imaginary. (2) It can be shown that the above expression for (<7ia[) 
reduces to (af) when T = 0, as it should. (3) The asymptotic behaviors of the summands 
when n — > oo, are ~ \ for A and ~ \ for B, hence both A and B converge. (4) If r 3> T, L, 
then A dominates over B, since the leading order of the summand in A is \fr while that in 
B is a constant independent of r as r — > oo. Thus (cici) ~ A. 

To proceed, let us now assume that r ^> T and r ^> L, then 



2rT + T 2 2rT , . 

p^V— ^-"v-p-* 1 (63) 

is a huge number. Thus for the third term in Eq. (|6~1~D, the sum should only start from n = p. 
We can now split the summation into two parts, i.e. terms with n < p and those with n > p. 
Using the asymptotic form of the summand for the part with n > p and defining m = [T/L], 
where [ ] denotes the integer part, one has 



/ <N~ r 2 ( 2L^ n , L*' 2rT-r- " 2r 2 \ 

Hence, it follows that 



,, ,.. r 2 /2L 4 * L l ' , 2rT-T 2 ~ 2r 2 \ 

1 ^ 1 < IiPZ» ( r 2 " S " + r 2 £ V " + + ? iH (65) 

Let us now evaluate the above expression term by term. One has, keeping in mind that 
p>l, that 



2rT-T 



2 V 



\l n2 + — 72 — ~ y n<2 +p 2 = pJ2y n2 /p 2 + 1 

n=l ' n=l n=l 



P 2 



f 1 Vx 2 + ldx^ -[v / 2 + ln(v / 2 + l)]p 2 , (66) 
j i/p 2 



and 



~ r 2 r 2 . . r 2 1 1 L L , 

V— = *(2,p) = r 2 . (67) 

^n 3 2 2 p 2 4 r T 1 J 



Here we have used Eq. (|63|) and the asymptotic expansion for ^(2,x) 
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where function \&(n, x) is defined as 

d n ^(x 



fy(n, x) 



dx n 



Noting that for p ^> 1, 

and letting 
we finally find 

k^;>i < 



1 



4- In r(x) 
ax 



71=1 Z 



ln(v^+l)], 



16tt 2 L 3 
Compare this result with 



L 4 
^^ 2 + 



if A 

2 p 2 



16tt 2 L 3 



2A + 



1\ LL 
2/7T 



IK) 



rC(3)r 2 



16tt 2 L 3 

It is seen that two successive pulses separated by T in time are uncorrelated (\(aia[ 
provided that 



IK 2 )«| 



or equivalently 



L 2 



rr, L 2 L T 

T> — = — L . 



However, if r <C L, one can show, by series expanding both A and B, that 

1 



16tt 2 L 3 



E 
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L 2 



)3/2 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



Clearly, in this case, |(<7ia{)| <C \ (&i)r\, if T > L, and |(<7i<7i)| ~ \(ct 2 )r\, when T < L. 

A few comments are now in order about the physical picture behind our correlation 
results. It is natural to expect from the configuration that the dominant contributions to 
the light cone fluctuation come from the graviton modes with wavelengths of the order of 
~ L. In other words, the lightcone fluctuates on a typical time scale of ~ 1/L. If the 
travel distance, r, is less than L, successive pulses are uncorrelated only when their time 
separation is greater than the typical fluctuation time scale. Otherwise they are correlated 
because the quantum gravitational vacuum fluctuations are not significant enough in the 
interval between the pulses. However, if r L 2 /T, then successive pulses are in general 
uncorrelated. Thus the correlation time for large r is of order L 2 /r, which is much smaller 
than the compactification scale L. We can understand this result as arising from the loss of 
correlation as the pulses propagate over an increasing distance. 
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4- Observational Limits 



Suppose that the experimental fractional resolution for a particular spectral line of period 
T is r. Then we must have 

^<r, (77) 

which, with Eq. (|5T)|), leads to a bound on L of 

L>^, (78) 

assuming L does not change over time. However, this bound can be trusted only when two 
successive wave crests are uncorrelated, when At is the expected variation in their arrival 
times. The most conservative constraint from this requirement is that L is smaller than the 
wavelength of the spectral line, T. Namely, 

^<T, (79) 
yielding a restriction on the range of spectral lines that we should use to get the lower bound 



Irtpi 



T>d-^. (80) 



If Eq. (BO) is approximately an equality, then Eq. (|78|) becomes 



i 




£>(^) 5 =m • (si) 

Obviously, the optimal lower bound would be deduced from the spectral lines of distant 
galaxies, possibly of cosmological distance, with the highest observed spectral resolution. 
For astrophysical sources of cosmological distance, spectral lines which satisfy Eq. ( |80| ) 
will have wavelengths > 1mm assuming a resolution of T w 10~ 3 and a cosmological travel 
distance. The detection of CO(l— ►O) line emission at 2.6 mm from luminous infrared galaxies 
and quasars provides the type of data needed to get a bound. According to Ref. |2~3 



the observed CO line resolution for the infrared quasar IRAS 07598+6508 , which is at a 
distance of 596 Mpc assuming H Q = 75 km s -1 Mpc -1 , is about T 10 -3 , leading to 

L>10 _1 mm. (82) 

Here, the size of the extra dimension has to be macroscopically large in order not to con- 
tradict the astrophysical observation. The lower bound given here is within the sensitivity 



of the recently proposed experiments for possible deviations from Newtonian gravity |p25 . 

Note, however, this bound could be pushed to an even higher value if we do not require 
L < T. This is in fact legitimate, since for astrophysical sources one usually has both r ^> L 
and r ^> T satisfied, and the condition for two successive pulses to be uncorrelated is really 

r>Y- ( 83 ) 
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Hence, we can use Eq. ([78]) to get a bound from experimental data as long as the resulting 
bound, L, obeys Eq. fl83"l), 

(84) 



L < VrT = L 



Thus, a stronger bound can be achieved, if we can find astrophysical sources of cosmolog- 
ical distance with observed spectra of much smaller wavelengths provided this condition is 
satisfied. The observation of 7 rays from astrophysical sources, such as gamma-ray bursters 
(GRBs) [H^L P rov ides such an opportunity. The use of these sources as probes of pos- 
sible quantum gravity effects has been explored by a number of authors p8|-p2| recently. 
Below we will select some of these 7 ray sources to calculate both the lower bound from 
Eq. flT8|) and L$ from Eq. (0). To be conservative, we shall assume a resolution of the order 
of unity for all the gamma rays we are going to consider, since we should at least have this 
resolution before we can talk with confidence about the observed energies (or frequencies) 
of the gamma rays. 

In the following Table, we list the source names, the observed frequencies u, the source 
distances D, and the calculated L& and L . 

Table: Bounds on L from GRB sources 



Source 


D(Mpc) 


v(Hz) 


Lb(mm) 


Lo(mm) 


GRB 930229 @ 


791 


4.8 x 10 19 


6.1 x 10 4 


1.2 x 10 10 


GRB 940217 


N 


385 


4.3 x 10 24 


2.6 x 10 9 


2.8 x 10 7 


GRB 930131 




260 


1.1 x 10 23 


4.6 x 10 7 


1.5 x 10 8 


Mrk 421 |HJ 


112 


4.8 x 10 26 


8.7 x 10 10 


1.4 x 10 6 


GRB 980703 


1592 


1.2 x 10 20 


3.1 x 10 5 


1.0 x 10 10 


GRB 980425 


Nl 


40 


5.4 x 10 20 


3.4 x 10 4 


8.3 x 10 8 


GRB 990123 |37| 


2400 


3.0 x 10 2U 


1.2 x 10 B 


8.5 x 10 9 



From this Table, one can find that the largest lower bound for L comes from GRB930131 
and GRB990123, which is 

L>10 7 mm. (85) 

For some sources, Mrk 421, for instance, one seems to get a much larger bound, which 
however can not be trusted, because the correlation condition Lf, < Lq is violated. The 
physical reason is that the frequencies for the gamma rays are so high that even travel over 
a cosmological distance does not wash out the correlation between successive wave crests. 

In principle, the above results only apply to the case where L is fixed. When L changes 
as the universe evolves, one should use Eq. (|53|), or Eq. fl55|), for the mean time deviation, 
At. In this case, we can set either a bound on L or a bound on the rate of change of L, if 
we can constrain either the rate of change of the size of extra dimensions over time or the 
present size L from other considerations. Let Lf be the size of the extra dimension at the 
present time and Lj be that at the time of primordial nucleosynthesis, and write 

^ = 1 + 5. (86) 
Li 
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According to Refs. pO| , pT[ , the observational limits, obtained by investigating the effects on 
the primordial nucleosynthesis of 4 He as a consequence of the time variation of fundamental 
constants such as the electroweak, strong, and gravitational coupling constants, imply that 
5 < 0.01. However, stronger limits on 5, which may be less reliable, arise from a detailed 
study of the events which took place 1.8 x 10 9 yr ago on the current site of an open-pit 



uranium mine at Oklo in the West African Republic of Gabon p2 |. This site gave rise to a 
natural nuclear reactor when it went critical for a period about 1.8 x 10 9 yr ago. The Oklo 



samples constrain the rate of change of extra spatial dimensions to satisfy the limits p3|j21 



< 1.9 x 10- 19 yr"\ (87) 



which translates to 5 < 10~ 10 . In either case, 8 <C 1. Thus all the results obtained so far 
for a fixed L hold for a changing L if renormalization with respect to L — > oo is adopted 
(refer to Eq. Q5BD). However if one chooses the renormalization with respect to the current 
size Lf, then combination of Eq. fl77|) and Eq. fl55|) gives rise to 

L > - >/5($) . (88) 

As a result, the following considerably smaller lower bounds can be deduced from Eq. ( |88D 
using the CO line data 

L>10" 2 mm, for 5 = 0.01, 

L>10~ 6 mm, for 5 = 10~ 10 , (89) 

and 

L>10 6 mm, for 5 = 0.01, 

L>10 2 mm, for 5 = KT 10 , (90) 

using GRB data. 

On the other hand, one can also use our result to place a bound on the rate of change of L 
if the present size of L can be fixed by other considerations. In this respect, it is well-known 
that the five- dimensional Kaluza-Klein theory provides an explanation of the quantization 
of electric charge, in the sense that all electric charges are multiples of the elementary charge 

e=— . (91) 
The corresponding fine structure constant is then 

«/ = % ■ ( 92 ) 

Setting aj to its present value, 1/137, we get an estimate of the size of the fifth dimension, 
L ~ 10~ 31 cm, in the original Kaluza-Klein model. For the case of renormalization with 
respect to L — >• oo, we can see from Eq. ([FT]), Eq. (|53"D and Eq. d86|) that we must have 
5 < 1, which is much weaker than the existing bounds on the change rate from primordial 
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nucleosynthesis and the Oklo samples. To discuss the case of renormalization with respect 
to Lf, the current size, let us write Eq. ( j88[) as 

Thus, using L ~ 10 -31 cm and the same CO line data as before, we find the following limit 
for change of the extra dimension in the Kaluza-Klein model 

5 < 1(T 58 , (94) 

or the rate of change by dividing 5 by the travel time, which is ~ 10 9 yr 



< 1(T 6 V • (95) 



This is much stronger even than the strongest bound arising from the observational limits 
on the time evolution of fundamental constants. 

To conclude, we have demonstrated, in the case of one extra dimension, that the large 
quantum lightcone fluctuations due to the compactification of the extra dimension require 
either the size of the extra dimension to be macroscopically large or rate of change of the 
extra dimension to be extremely small. This result seems to rule out the five dimensional 
Kaluza-Klein theory, or at the very least, place strong limits on the rate of change of the extra 
dimension. We must point out that the rate of growth of At with r depends crucially on the 
number of spatial dimensions. In four dimensions, At oc y/r, while in five dimensions At oc r. 
One expects that in larger number of dimensions, there will be an effect of compactification, 
but its details need to be determined by explicit calculations for particular models. This is 
the topic for the next section. 



B. Higher dimensional models 

There was no real reason to extend the Kaluza-Klein idea beyond five dimensions until 
the emergence of non-abelian gauge field theories which have had a profound impact on 
theoretical physics since their invention by Yang and Mills in 1954. It was suggested by 
DeWitt [^U as early as in 1963 that a unification of Yang- Mills theories and gravitation could 
be achieved in a higher dimensional Kaluza-Klein framework. Nowadays, the possibility of 
unifying all the known interactions in nature in higher dimensional spacetimes has been 
actively pursued in the context of eleven dimensional supergravity and ten dimensional 
superstring theories. A necessary ingredient of all these higher dimensional models is the 
compactification of the extra dimensions to a very small size so as to leave the ordinary 
four- dimensional "large" world. 

In this section, we examine, case by case, the lightcone fluctuations arising from quantum 
gravitational vacuum fluctuations induced by the periodic compactification of extra, flat 
spatial dimensions in higher dimensional models up through 11 dimensions and make a 
conjecture about arbitrary dimensions. 
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a. The case with n = 2 This is the 6 dimensional spacetime with 2 extra dimensions. 
We obtain, after performing the differentiation in both F xx and H xxxx (See Eqs (|37|) and 



!))• 



with 



G xxxx (t, x, 0, t', x', 0, 0, mi-L, m 2 L) = A 2 (p) + B 2 {p) 



A 2 ( P ) 



471-3 



54p 6 



2^4 



+ 



102a^p 



4^2 



3a 2 p- 



and 



£ 2 (p) 



47T 3 
hi 



6p 7 



(p 2 + ai 



(p 2 + a 2 2 ) 5 (p 2 + «|) 5 2(p 2 + ai) 5 



27a 2 p 5 27a^p 3 

~ ' Tl a ON 11 I" 



(96) 



(97) 



3a 2 p 



(p 2 + Oi\ 



4(p 2 + a|)t 8(p 2 + «I) 



yV + Qj + p " 
s /p 2 + a 2 2 -p y 

where we have introduced a parameter 



(98) 



a, 



^ 2 E 



m 2 



(99) 



i=l 



Upon noting that G xxxx is an even function of p, it follows that, 



rb rb 2 +oo „ r 

/ dx dx'G^ xx (t,x,0,t , ,x , ,0) = ]J 2(r-p)(A 2 + B 2 )dp. 

Ja Ja i=l mi =-oo J ° 



(100) 



Performing the integration (integrate by parts for those terms involving logarithmic func- 
tions), we find 



2 +oo „ r 

1] E / 2(r-p)(A 2 + B 2 )dp 

i=l m;=-oo J0 

7 2 (m 2 + m 2 ,) 



1 

2vr 3 L 2 



oo ' oo ' 

E E 

mi=l m2=l 



+ 



10 7 4 



(7 2 + m 2 + m 2 ) 3 3(7 2 + m 2 + m 2 ) 3 
87 s 8 7 5 + 4(m x + m 2 )7 3 + (m\ + m 2 ) 2 7 



3(m 2 + m 2 )(7 2 + m\ + m 2 ) 3 



2(m? + mi + 7 2 ) 7 /2 



x In 



mf + m 2 + 7 2 + 7 
m 2 + m 2 + 7 2 — 7 / 



101) 



where 7 = r/L is a dimensionless parameter. The above double summation is by no means 
easy to evaluate. However, we are interested in the case in which the travel distance r is 
much greater than the size of the extra dimensions L, i.e., 7 ^> 1. It then follows that the 
summation is dominated, to the leading order, by 
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2tt 3 L 2 



E E 

mi=l m2= 



87 s 

1 3(m 2 + m^)^ 2 + m 2 + m 2 ) 3 ' 



(102) 



which can be approximated by the following integral when 7 1 



3tt 3 L 2 hh 



roo roo 

/ dxi \ 

7i/ 7 A/ 7 



1 



(a; 2 + a;|)(a;f + x| + l) 3 3tt 3 L : 



111(7) > as 7 ^ 00. (103) 



An easy way to see the above behavior is to note that the contribution to the integral is 
dominated by the region around (1/7,1/7) since the integrand dies away very quickly as X\ 
and xi increase, and then change to polar coordinates to evaluate the integral while using 
the fact that around (I/7, 1/7) the integrand is approximated by 1/ (x\ + x\). Therefore one 
finds that 



R 



Hi) , 



(104) 



6tt 3 L 2 

and the mean deviation from the classical propagation time due to the lightcone fluctuations 



At = 



l327vG e 
6tt 3 L 2 



At 



P i 



3tt 2 



71 



ln(r/L) . 



(105) 



Here we have used the fact that G± +n = G^L n . Note that here the At dependence on 7 is 
square root of a logarithmic function which is quite different from linear dependence in the 
case with n = 5 discussed above. So, the mean time deviation here grows much more slowly 
as the travel distance r increases or the compactification scale decreases. As we shall see 
later, this seems to a general feature for n > 5. 

b. The case with n=3 This is the 7 dimensional spacetime with 3 extra dimensions. 
The relevant two-point function can be shown to be 



G xxxx (t, x, 0, t', x', 0, 0, miL, ...,m 3 L) 



8tt 3 



24p 8 



336a3p 6 



3^4 



+ 



280« 3 p' 



« 3 (P 2 + «|) 6 (p 2 + «|) 6 2(p 2 + a 3 )6 



(106) 



which leads to 
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/ dx I dx G xxxx (t, x, 0, t , x , 
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3 00 
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m, s 1 ^3(EL 1 m 2 )f( 7 2 + ELi^) 4 
56 7 6 



3(El 1 m 2 )n7 2 + EL^) 4 



(107) 



The triple summation is dominated, to the leading order, by the first term when 7 ^> 1, 
which again can be approximated by integration to be 



3"7T 3 L 3 Jl/7 



dxi 
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1/7 
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(108) 
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as 7 — > oo. Thus the mean deviation from the classical propagation time due to the lightcone 
fluctuations 



At 



R 



' 32nG 7 
12tt 3 L 3 



ln( 7 ) 



3tt 2 



ln( 7 ) 



>/Mt)- 



(109) 



c. The case with n=4 This is the 8 dimensional spacetime with 4 extra dimensions. 
The relevant two-point function is 

G xxxx (t, x, 0, t', x', 0, 0, miL, m A L) = A 4 (p) + B A (p) 



(110) 



with 
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(111) 



and 
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8vr 4 
In 



240p 7 



4^3 



540a|p 5 90a|p 

13~ + 



+ 



Yba\p 



(p 2 + al)- (p 2 + al)- (p 2 + a\)- 4(p 2 + al)~ 
p 2 + aj + p 



x 



^p 2 + aj -p, 
One finds after carrying out the integration 
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The summation is seen to be dominated by the second term when 7 ^> 1 and in that case 
the summation turns out to be approximated by an integral as 
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Therefore, one obtains 
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d. The case with n=5 This is the 9 dimensional spacetime with 5 extra dimensions. 
One finds in this case that 
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and 
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The summation, to the leading order, can be approximated by integration 
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as 7 — > 00. Thus the mean deviation from the classical propagation time due to the lightcone 
fluctuations 
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e. The case with n=6 This is the 10 dimensional spacetime with 4 extra dimensions as 
motivated by superstring theory The relevant two-point function is given by 
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One finds after carrying out the integration 

rb r-b 

\ dx dx'G* xxx (t,x,0,0, lf,a/,0,miL) 

J a J a 

424 7 6 



1 



6 oo ' 



IIE 



687 s 



3907 4 



+ 



195 7 2 (EJLi m i 



4 7 



10 



(7 2 + ELi (7 2 + Eti m?) 6 3(ELi ™?) 3 (7 2 + E?=i ™ 2 ) 5 

400 7 5 + 100(Eti "i^T 3 + 15(Eti m?) 2 7 



2(Eti ™ 2 + 7 2 ) 11/2 



In 



E-=i < + 7 2 + 7 



(124) 



The summation is dominated by the second term when 7 ^> 1 and thus is approximated by 
an integral as 
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/. The case with n=7 This is the 11 dimensional spacetime with 7 extra dimensions . 
The relevant two-point function can be shown to be 
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Thus we find that in all of these models in which there is more than one flat compactified 
extra dimension, the lightcone fluctuation effect grows only logarithmically with distance. 
Hence we are unable to derive any constraints on these models. We have not been able to 
prove that this behavior holds for any number of flat extra dimensions greater than one, but 
conjecture that this is the case. 



IV. PARALLEL BRANE- WORLDS SCENARIO 

In this section, we examine the lightcone fluctuations due to the presence of two 3+1 
dimensional hyperplane boundaries, i.e., " 3-branes", living in extra dimensions separated 
from each other by some distance. This framework is motivated by the recent proposal to 
resolve the unnatural hierarchy between the weak and Planck scales 0. In this scenario, 
four- dimensional particle theory, such as the Standard Model, is confined to live in one of 
the branes, but gravity is free to propagate in the higher dimensional bulk. Therefore the 
bulk spacetime is dynamical, and the 3-branes can not be rigid, but must undergo quantum 
fluctuations in their positions, which we assume to be order of l p , the Planck length in higher 
dimensions. To be more specific, we suppose one brane is located at the origin and the other 
at (0, 0, 0, zi, Z2, z n ), and we shall examine the effect of lightcone fluctuations due to 
gravitons in the bulk by looking at a light ray traveling parallel to one of the boundaries, in 
the x— axis, for example, but separated from it by a distance z ~ l p . This feature is to reflect 



the quantum fluctuations in the position of the brane | 55fl . We consider both the Neumann 
and Dirichlet boundary conditions for the graviton field at each brane hyperplane. Once 
we have the graviton two-point functions without any boundary, those with two parallel 
plane boundaries can also be found by the method of image sum. If gravitons satisfy the 
Dirichlet boundary condition, then the renormalized graviton two-point function is given by 
the following multiple image sum 

Gfj, vpa (t, Z ii t > Z i) 



n +00 
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Again the prime denotes omitting the m; = term in the summation. For the Neumann 
boundary, the renormalized graviton two-point function becomes 

G pup(T (t, Zi, t , Zj) 
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= II J2 G^upait, z h t', 2miL + z[) + G^patt, z i7 t', 2miL - z[) 

i=l rrii=l 

+G IJ , upa (t, Zi, t', -2miL + z'i) + G pvpa (t, z h t', -Im^ - z[) 

+G pupa (t, Zi, t' , —z'j) . (132) 

Let us examine a light ray propagating along the x-axis starting from (a, 0, 0, zi, z n ) 
to (b, 0, 0, zi, z n ). For simplicity, let us assume that z\ — ■ ■ • — z n — z ~ l p j y/n. Then 
the mean squared fluctuation in the geodesic interval function is 

K 2 ) = \{b-af fdx fdx' G^ xx (t,x,0,0,Zi, if, x', 0, 0, Zi ). (133) 

O J a J a 



A. Five dimensional theory 

Here we have one extra dimension, and the relevant two-point function is 
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where Azi = Z\ — z[. We then have 
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If gravitons satisfy the Dirichlet boundary conditions, it follows that 
K 2 ) = \r 2 fdx fdx' G^ xx (t,x,0,0, Zl , t',x',0, Zi ) 

o J a J a 

1 1 oo 

= ~-r 2 f(2z, r) + -r 2 £ [2/(2mL, r) - f(2z - 2mL, r) - f(2z + 2mL, r)] . 

8 8 m=l 

(136) 

Here we are interested in the case in which r z and L~^> z. One then finds that 

* -132^7 = -l3^(f) 3 ' (137) 

and 
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Here £(5) is the Riemann-zeta function. Consequently, we find 
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thus, 
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If instead of the Dirichlet boundary condition, gravitons are forced to satisfy the Neu- 
mann boundary condition, one finds 
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Here we find that the choice of different boundary conditions has very little effect on the 
growth of the lightcone fluctuations as long as L ^> z, although the fluctuations are slightly 
larger in the case of the Neumann boundary condition. Comparing this At with that in 



the case of periodical compactification Eq. (|50[) , one can see that the effect of light cone 
fluctuation here is larger for a given size of extra dimensions, L. Another very distinctive 
feature is that here At increases as L increases provided that L and z are independent of 
each other, which is in sharp contrast with the case of periodical compactification. As before, 
the fluctuation in the flight time of pulses, At, can be applied to the successive wave crests 
of a plane wave. This leads to a broadening of spectral lines from luminous sources. In case 
of periodical compactification, we used the data from gamma ray bursters to derive a very 
strong lower bound on L. In the present framework, however, we have an intrinsic lower 
bound L > z. Recall that z is the Planck length in higher dimensions which may be much 
larger than that in 4 dimensions. In the picture of Ref. pfl, z and L are not independent, 
and may be expressed in terms of the Planck mass in 4 + n dimensions, M^ +n ^ , as 
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Let M^ 4+ra ) be ~ ITev. Then for an astronomical source of cosmological distance, such as 
gamma-ray bursters with redshift ~ 1, or r ~ 10 28 cm, we get the following estimate for the 
n — 1 model 



At ~ 10 60 x l(T 44 s = 10 16 s. 
That is far too large, so that the n — 1 model can be ruled out completely. 
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B. Higher dimensional theories 

We shall only study the 6 dimensional case in some detail, but shall list results for other 
cases. It is not difficult to see that the relevant two-point function can be obtained from 
Eq. ( p6|) by replacing a 2 there with J2i=i ^ z f- Therefore one easily finds from Eq. ( |101| ) 
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Then it follows that 



= g r2 dx dx' G xxxx (t, x, 0, 0, z, z, t 1 ', x', 0, z, z) 
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Here we are interested in the case when r ^> z and L ^> z. One finds that 
and 
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If we assume the travel distance r is much larger than the size of the extra dimensions L, 
i.e 7 > 1, g(r, z, L) can be approximated by integrals as follows 
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Therefore, 
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Hence, 
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Notice here that the leading term does not depend on the travel distance r and moreover 
the next order corrections decrease as r increases. Recall that 7 = r/L. Interestingly, the 
lightcone fluctuation grows as the compactification scale increases. 

However, if we change to the Neumann boundary condition, the behavior will be different 
as we can see from the following analysis. Let us note that in this case 
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Notice the sign changes in the above expression as compared to Eq (|150|) . Then it follows 
that 
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and 
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Here we have two different contributing terms to the mean deviation in time; one is inde- 
pendent of r but increases as L increases, the other grows logarithmically as r increase or L 
decreases. 

Similarly, we have calculated cases up to n = 7 as motivated by string/M theory. The 
results are (all in units of t p i, the Planck time in four dimensions): 
For n=3, 
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For n=4, 
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For n=5, 
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For n=6, 
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For n=7, 
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A few comments are now in order for all the cases we have examined: First, to the 
leading order At does not grow as r increases for the Dirichlet boundary condition. Second, 
At behaves differently for different boundary conditions, but the logarithmic dependence 
associated with the Neumann boundary condition can be neglected as far as observation 
is concerned because logarithmic growth is extremely small. Third, our results seem to 
suggest a leading order behavior of (L/z) n ^ 2 for any dimensions. Finally, our results are 
quite different from those obtained by Ref. |^0[ where only the contribution of the F term 
was considered and only one dominant term in the sums computed. Our results are much 
smaller than those obtained by these authors because of cancellations among the various 
terms. 

We wish to gain an understanding of how large the effect of lightcone fluctuations can 
be in these higher dimensions and whether they might be observable. 
From Eq. (|145|) and the above results, one finds 
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This result reveals that the smaller is M^ +n ^> , the larger is the effect of lightcone fluctuations. 
According to Ref. 0, M^ +n) may as low as the order of one Tev. This leads to 

At~l(T 2 V (165) 

This is a tiny effect by conventional standards. However if we note that the above mean 
time deviation is equivalent to an uncertainty in position 

Ax~l(T 21 m (166) 

and the operation of gravity-wave interferometers is based upon the detection of minute 
changes in the positions of some test masses (relative to the position of a beam splitter), we 
can see that this effect might be testable in the next LIGO/ VIRGO generation of gravity- 
wave interferometers pETP^fl. It constitutes an additional source of noise due to quantum 



gravity. Currently, the sensitivity of these gravity-wave interferometers has already reached 
the order of lCT 19 m ||43|| . Note that Amelino-Camelia [H and Ng and van Dam [[45] have 



also proposed rather different quantum gravity effects which might also be detectable by 
laser interferometers. 

A more complete discussion of the observability of the At given by Eq. (16~3) should 



involve a calculation of the correlation time, analogous to that performed in Sect. |III A3| 
for the five dimensional Kaluza-Klein model. This calculation has not yet been performed. 
However, it is reasonable to guess that the result will be of the order of or less than Tev scale 
which characterizes this model. Recall that in the five dimensional compactified model, the 
correlation time was found to be much smaller than the compactification scale when the 
travel distance is large. If this guess is correct, then the correlation time is much smaller 
than At itself. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have examined the effects of compactified extra dimensions upon the 
propagation of light in the uncompactified dimensions. There are nontrivial effects that 
arise from quantum fluctuations of the gravitational field, induced by the compactification. 
These effects take the form of lightcone fluctuations, variations in the flight times of pulses 
between a source and a detector. The crucial quantity describing these fluctuations is At, 
the expected variation in arrival times of two successive pulses which are separated by more 
than a correlation time. In Sect. [11 B|, we gave a derivation of the formula for At using 



the geodesic deviation equation. This derivation allowed us to discuss issues of gauge and 
Lorentz invariance. In particular, it demonstrates that At is gauge invariant. All of the 
subsequent explicit calculations are performed in the transverse-tracefree gauge. 

As a prelude, we found the graviton two-point function in this gauge in Minkowski 
spacetime of arbitrary dimension. The two-point function in a higher dimensional flat com- 
pactified spacetime is given as an image sum. We then calculated At in the five-dimensional 
Kaluza-Klein model, five-dimensional flat spacetime with one periodic spatial dimension. We 
found that At grows linearly with increasing distance between the source and the detector, 
and is inversely proportional to the compactification length, L. This result differs from the 
square root dependence on distance that was found in four-dimensional flat spacetime with 
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one periodic spatial dimension ||11|| . This demonstrates that lightcone fluctuation effects are 
rather model-dependent. We also calculated the correlation time in the five-dimensional 
model and found that it is typically small compared to the compactification scale L. This 
allows us to place very tight constraints on the parameters of this model. 

We favor the viewpoint that the lightcone fluctuation effects should vanish only in the 
limit that L — > oo. If one adopts this view, then the five- dimensional Kaluza-Klein model 
can be ruled out. Data from gamma ray burst sources imply a lower bound on L which 
is larger than upper bounds obtained from other considerations. However, another logical 
possibility is that lightcone fluctuation effects happen to vanish ((erf) = 0) at the present 
compactification scale. Even if one adopts this viewpoint, one still obtains very strong 
constraints on the fractional change in L which can have occurred over a cosmological time 
scale. This in turn tightly constrains any five-dimensional Kaluza-Klein cosmology. 

We also examined an alternative five- dimensional model (the brane worlds scenario) in 
which gravitons satisfy Dirichlet or Neumann boundary conditions on a pair of parallel four- 
dimensional branes (one of which represents our world). In this model, we again find that 
the lightcone fluctuations are so large that the model can be ruled out. 

We next turned our attention to models with more than one extra dimension. In the 
case of two or more flat, periodically compactified dimensions, we found that At grows only 
logarithmically with distance, and that no constraints may be placed on these models. In 
the case of the brane worlds scenario with more than one extra dimension, we found that 
At approaches a constant which can be of the order of 10~ 29 s. This would produce a source 
of noise in laser interferometer detectors of gravity waves, which may eventually be within 
their range of sensitivity. 

In summary, compactified extra dimensions have the possibility to produce observable 
effects by enhancing the quantum fluctuations of the gravitational field. These effects might 
be used either to place constraints on theories with extra dimensions, or else possibly even- 
tually to provide positive evidence for the existence of extra dimensions. Although our 
results for more than one flat compactified extra dimension are too weak for either of these 
purposes, the required calculations for models with curved extra dimensions have not yet 
been performed. Another model which has not yet been examined in this context is the 



Randall- Sundrum model ||46|| , with an uncompactified fifth dimension. In this latter model, 
propagating graviton modes are effectively confined within a finite volume in the fifth di- 
mension, so one might expect nonzero lightcone fluctuations. 
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VI. GRAVITON TWO-POINT FUNCTIONS IN SPACETIME WITH 
ARBITRARY NUMBER OF EXTRA DIMENSIONS 



Here we evaluate the functions D n (x.x'), F™Ax,x') and H^ kl (x,x') defined in Eqs. 



(j2q) and (f28|) , respectively. Once these functions are given, the graviton two point functions 
are easy to obtain. Define 

i?=|x-x / |, At = t-t', k=\k\ = w, (Al) 

and assume n extra dimensions, then 

D n (x, x') = [ ^!!!^ e ^-(x-x') e — (t-f) 

Tip poo pir 



poo pit 

/ k n+1 e- ikAt dk / d0 1 sin 1+n ^e 
Jo Jo 



2(2vr) 3 + n Jo Jo 

pTT p7V p 

x / d0 2 sm n 9 2 ... I d6 n+1 sin6» n+ i / 
Jo Jo Jo 



2tt 

dd, 



n+2 



poo pi 

J k n+1 e~ ikAt dk J e ikRx {l - x 2 ) n/2 dx 



a n Re f°° 
2(2vr) 3+n 

■p^- r k n+1 e- ikAt dk [\l-x 2 ) n / 2 cos(kRx)dx 
(27r) 3+n Jo Jo 

a n ^2^r(| + 1) Re 1 



2(27r) 3 +™ Jo 



poo i i 

/ k^Jn+i(kR)e- ikAt dk 
Jo 2 



a n V7r2^r(f + 1) Re 1 . r°° n+i 



poo , - 

lim / k—J n±1 (kR)e- ak dk 

0++iAtJo 2 



\3+n D «±l 



2(27r) 3 +™ a^0++iAt Jo 

a n 2«r(f + 1) 2 i r(f + 1) i 



(2tt) 3 +™ (i? 2 - At 2 )™/^ 1 4tt^ (i? 2 - At 2 )™/^ 1 



Here we have defined 



i"K Pir r2ir 2,71^^^ 

i n = de 2 sm n e 2 ... d6 n+1 sm6 n+1 d6 n+2 = —- (A3) 

JO JO JO 1 n + 1 



r(f + i) 



and used 



/ cos(kRx)(l - x 2 ) n/2 dx = y^T(- + 1) ( -=- ) 2 J^+i(kR) , (A4) 
Jo 2 2 \kRJ 2 



and 



(A2) 



e~ ax JJf3x)x v dx = (WEC^+V^ to > -1/2. (A5) 
V 7 v^F(a 2 + /5 2 )^ 7 V 7 

When n is odd, D n (x,x') should be taken to be zero when R 2 < At 2 . 

Let us now turn our attention to the calculation of i 7 ^- and HjM. We find 
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F n( >\ _ Re f j3+n k ^7 ik-(x-x') -iw(t-f) 

^^' XJ ~ (2vr) 3 W 2cu 3 

= —^—di& / k^e'^dk / d0! sin 1+ "# ie ^ cosei 
2(2tt)3+™ jy Jo 

/■7T /"7T /*27T 

x / sin n 6* 2 ... / d# n +i sin0„ + i / dO n+2 
Jo Jo Jo 

= anRe d t d' f°° k n - l e- lkM dk I'' (1 - x 2 )"/ 2 cos(kRx) dx 
{27l)- i+n J Jo Jo 



^d' f — \+r [ k^Jn+i (kR)e 
J \R— Jo 2 



-ikAt 



3+n 



-Re , ( n — 1 1 



^(^^r r k^J^{kR)e~ ikAt dk 
J \ R 2 R— Jo 2 



^ . v 3+n * 7 I 7-)9 „n— 1 / 

2(2tt)— j V i? 2 i?— Jo 

' f ° k^ Ju^(kR)e- ikAt dk ) , (A6) 



2 



where we have utilized a recursive formula for Bessel functions 

zJ v ^{z) + zJ^+i^) = 2vJ v {z) . (A7) 

Similarly, one finds that 

To proceed further with the calculation, we need to deal with the cases when n is odd or 
even separately. 



A. The case of odd n 



Assume n = 2m + 1 and define 

Re r°° 



l-fp poo 

S ^ = R^J k m - 1 J m+1 (kR)e- lkAt dk, m>0, (A9) 



Tip roc 

T ( m ~ !) = I k m - 1 J m ^ 1 (kR)e- lkAt dk 
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Re 



R r 



- lim / k 171 - 1 J m JkR)e- 

1 a^O++iAt Jo 



■alphak 



dk 



2 m " 1 r(m - 1 /2) VR 2 - At 2 
7^ # 2 (# 2 - At 2 ) 

(2m -1)!! Vi? 2 - At 2 
(2m- 1) i? 2 (i? 2 - At 2 )" 1 ' 



m > 1 



where we have appealed to integral (6.623.1) in Ref. |47|.The above result holds for R 2 
and T(m — 1) is zero when R 2 < At 2 . Then it follows from Eq ( |A6| ) that 



and 



2/77, 

ST^) = ~iwS(. m ~ i) — — i) 
Using the recursive relation Eq (|A12|) , we can show that 

b i m ) ~ R 2m ( 2Jfe)!! i? 2m - 2fc 



^(0) 
(2m)!! 



(2jfe-l)!! i? : 



2A- 



i - V 

^ (2fc)!!(2fc- 1) (i? 2 - At 2 ) fc 
(2Jfe + l)!! R 2k 



# 2m (2Jfe)!!(2Jfe + l)(2Jfe - 1) (R 2 - At 2 ) k ' 



Here 



5(0) 



^3 for/? 2 > At 2 
fori? 2 < At 2 



as reported in Ref. ||12|| . For the sake of completeness, we give its derivation below 

J^kR) cos(fcAt) 



S(0) = % lim 

R a^0++iAt 



k R Jo 



k 



i cos(arcsin(At/i?)) = Vfl2 ~ At2 fori? 2 > At 2 



fori? 2 < At 2 . 



dk 



If we define 



Q(m) 



Re 



R r 



- / k m - 3 J m+1 (kR)e~ ikAt dk, m>0 
1 io 



then it is easy to see that 



Tj2m+1 
n ijkl 



2(2n] 



m+2 1 3 



d i d , j d k d , l [Q{m 



39 



and 



2m 1 
Q(m) = —Q(m-l)-—S(m-2) 



(A18) 



The above equation applies for m > 2. To use it to get a general expression, we need Q(0), 
which is given by 



1 /» 1 

Q(0) = - / j^JiikR) cos(kAt) dk 
R Jo k 6 



lim 



Jt(kR) cos(kAt) dk 



£-o R Jo (k 2 + (3 2 ) 

1 e-P At h(f3R) 1 1 



lim 

P-+o R 



ff 1 



2(3 2At 



(A19) 



This leads to a vanishing H^b We next need Q(l), which can be calculated, using integral 
(6.693.5) in Ref. |47j, as follows 



1 [™ J 2 (Rk) cos(Atk) 



R 2 



o 



k 2 



dk 



i 

R 2 





cos(arcsin(At/i?)) + § cos(3 arcsin(At/i?))l for R 2 > At 2 , 

fori? 2 < At 2 . 



1 _ At 2 \ VR 2 -At 2 
3 3R 2 J iJ 2 







(j-gf)S(O) fori? 2 > At 2 
fori? 2 < At 2 



In the above calculation, we have made use of the following trigonometric relations 

cos(3a;) = 4cos 3 (x) — 3cos(x), 



cos arcsmx 



vT 



(A20) 



(A21) 



Therefore one finds, using the recursive relation Eq ( [A18| ). 
WK) 2R 2 ™~ 2 ^ { ) i? 2 , 2 ^ (2k)\\ R 2m - 2k 



(2m)!! 
R 2 



m-2 



\Q{i) 



(2j + 1)!! 



i? 2 ^ 



m k—2 

+ 21(2^- 2)(2j)!!(2j + l)(2j - 1) (i? 2 - At 2 > 



r 5(0) 



This expression can be simplified if we note that 

1 ^ 1 



fc=2 !y 



E 

fc=2 



1 m — j — 1 

k(k - 1) ~~ m(j + 1) 



(A22) 



(A23) 
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and 



m k—2 

EE 

k=2 j=0 

m—2 m 

= E E 



(2j + l)\\ 



R 2 i 



5(0) 



2k(2k - 2)(2j)!!(2j + l)(2j - 1) (i? - At 2 )-? 
(2j + 1)!! 

j=0 k=j+2 2k(2k - wmm + Wo - 1) (r 2 - At 2 ) 

(m-j-l)(2j + l)H 



7 5(0) 



E 



Am(j + 1)(2j)!!(2j + l)(2j - 1) (i? 2 - At 2 ) 



5(0) 



So, we have in this case 



(A24) 



D 



2m+l 



(2m+l)!! 1 



fori? 2 > At 2 , 
fori? 2 < At 2 , 



(A25) 



Tp2m+l 



i oo> fMllc^y ( 2fc + 1 ) !! fi2fc ' 

2 ( 2vr )m+2 » j ^ ffm ^ ( 2 jfc)!!(2A; + l)(2Jfe - 1) (i? 2 - At 2 ) fc , 



, (A26) 



and 



ZT2m+l 
n ijkl 



2(2tt 



\m+2 



,f(2m)!! 



i? 2 



m—2 



m-2 



+ E 



(m-j -l)(2j + l)! 



i? 2 ^ 



4m(j + l)(2j)!!(2j + l)(2j - 1) (i? 2 - At 2 > 



-5(0) 



(A27) 



B. The case of even n 



Let n = 2m with m = 1, 2, 3.... The graviton two-point functions for m = correspond- 
ing to the usual 4 dimensional spacetime have been given previously [|ll|]. The analog of 
Eq. ( [A 11| ) for this case is 



Here 



F?™ = r di& (S(m - -) . (A28) 

3 2(27r) m+ l 3 \ 2') v ' 



S(m - 1/2) = -^^(m - \) ~ T(m - |) . (A29) 

Using this recursive relation, we can express S(m — 1/2) in terms of S(l/2) which is calcu- 
lated, by employing 

(ABO 
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to be 



S(l/2)- 




k~ 1/2 Js(Rk)cos(Atk) dk 
1 r°° d fsm(Rk)\ , . , s „ 



sin(Rk) cos(Atk) 



k 



At 
4R3 



In 



At 
+ ~R? 



sxn(Rk) sin(Atk) 
k 



dk 



It then follows that 

. . (2m - 1 
S(m- 1/2) = 



R + At 



R- At J 



(2m-l)\\T(k 



(A31) 




T(k-V 



R 



2k 



(R 2 



Similarly, one has for H, 



2m 
ijkl 



H 



2m 
ijkl 



2(2tt)' 



and 



9rn _ 1 3 1 5 

Q(m - 1/2) = ~^Q{m --)- — 2 S{m - -) 



At 2 )*- 1 
(A32) 

(A33) 
(A34) 



Now the calculation becomes a little tricky. First, let us note that H^ kl has already been 
given [IT | and the recursive relation Eq ( |A34j ) can only be applied when m > 3. So, we 
need both H 2 - kl and Hf- kl or Q{l/2) and Q(3/2) as our basis to use the recursive relation 
for a general expression. Because there is an infrared divergence in the Q(l/2) integral, so, 
as we did in the 4 dimensional case, we will introduce a regulator (3 in the denominator of 
the integrand and then let (3 approach after the integration is performed. Noting that 



Ja(z) 

2 



smz 




cos z 



(A35) 



we obtain 



0(1/2) = i 
Rz 




k 2 Js(kR) cos(kt) dk 



dk , „ . 1 
— sin kR cos kAt — — - 
k R l 



dk 

— cos kR cos kAt 

k 6 



+ 



1 1 d 
R^^dpJo 



d r°° sin kR cos kAt 

2(3 d(3 Jo k 2 + (3 2 

00 k cos kR cos kAt „ 
- dk 



dk 



k 2 + (3 2 



(A36) 



42 



We next use 

I 132 + J dx = ^e-^{e b ^Ez[(3(a - &)] + e^Ei^d + b)}} 

-j3^{e^Ei[-P{a + b)} + e~ bf3 Et[- (3{a - b)]} , (A37) 

l X XCOS if )C ° 2 S(M dx = ~e-^{e b ^Ei[(3(a - b)} + e'^Ei^a + b)}} 
Jo p -\- x 4 

-\ e ^{ e ^Ei[-P{a + b)} + e- hp Ei[-(3{a - b)}} , (A38) 

where Ei(x) is the exponential-integral function, and the fact that, when x is small, 

Ei(x) « 7 + In |x| + x + ^x 2 + -^x 3 + 0(x 4 ) , (A39) 

where 7 is the Euler constant. After expanding Q(l/2) around (3 — to the order of /3 2 , 
one finds 

At 2 At /At 2 \ /R + At\ 2 \ . . 

Note, however, that what we need is Hijki which involves differentiation of (5(1/2), therefore 
we can discard the constant and divergent terms in (5(1/2) as far as H^ki is concerned. To 
calculate (5(3/2), let us recall that 

Q(3/2) = ^Q(l/2) - ±S{-l/2) (A41) 



and note that S(— 1/2) is given by y2/-K times Eq (A19) in Ref. |Tl] Thus, we have 

^/ in /2/ 1 At 2 At /At 2 \ /i? + At\ 2 \ /A . 

With (5(3/2) at hand, it is easy to show that for an arbitrary m > 3 

Q(m - 1/2) - (2m ~ 1)!! C(3/2) -±Y ( 2 ^- 1 ) !! ^- 5 / 2 ) 
Wlm L/2) - 3B?™~ * Q{6/ Z) JPj^ (2*-l)!! i? 2m " 2fc 

^ " ( ^(3/2) - E _ -4—5(1/2) 



R 2m-4 V3^ W y ^ 3 (2fc- l)(2fc-3) 
1 /2~™^ 2 2 1 '- 2 r(j-l) i? 21 ' 



i? 4 V vr ^3^2 (2j - 1)!!(2A; - 1)(2* - 3) (i? 2 - Ai 2 )^ 4 
(2m - 1)!! / 1 



R 2m-A V 3^( 3 / 2 ) E (2fc-l)(2A;-3) 5(1/2) 

1 Fl^ (m — j — l)2^- 2 r(j - 1) R 2 i 
+ B iS i^p 2 (2m-l)(2j + l)!! {R 2 - At 2 )i- 1 



(A43) 
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Here in the last step, we have made use of the following results 



1 



1 



j'+i 



1 



Jr +2 (2* - 1)(2* - 3) § ( 2 * - W 2 * - 3 ) £ (2* " 1)(2* - 3) 

m — j — 1 



(2m -1) (2J + 1)' 



(A44) 



and 



m, fc— 2 m k—2 m—2 m 

E E /OX*) = E E /OX*) = E fU) E • 

fc=3j=2 fc=4j=2 i=2 fc=2+j 



(A45) 



Consequently, we obtain 



D 



2m 



2 m m\ 



1 



(27r) m + 2 (i? 2 - At 2 )"^ 1 ' 



(A46) 



2m 



_J_ f (2m-l)H 

(27r) m + 2 4 H R 2m 



AR \R-At) 



R 



2k 



1 ™ 2 fc " 2 r(A; - 1) 
"i? 2 ^ (2fc-l)!! (if 2 -At 2 )*- 1 



k=2 



(A47) 



and 



Tj2m 
n ijkl 



1 



(27r) m + 2 ' ' 
At 2 



At /At 2 
24R3 V"^ 2 " 



- 1 In 



(R + At\ 2 
\-R-At) 



18R 2 



6i? 4 



-E 



fc=3 (2fc-l)(2fc-3) 
1 ^ 2 (m - j - l)2^'- 2 r(j - 1) R?i 



i? 2 4i? 3 n v J R-At; 
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